ON THE BOUNDARY OF THE PSEUDOSPECTRUM 
AND ITS FAULT POINTS 
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Abstract. The study of pseudo-spectra of linear transformations has become 
a significant part of numerical linear algebra and related areas. A large body 
of research activity has focused on how to compute these sets for a given 
spectral problem with, possibly, certain underlying structure. The theme of 
this paper was motivated by the question: How effective are path-following 
procedures for tracing the pseudospectral boundary? The present study of the 
mathematical properties of the boundary of the pseudospectrum is the result. 
Although this boundary is generally made up smooth curves, it is shown how 
the Schur triangular form of the matrix can be used to analyse the singular 
points of the boundary. 

1. Preliminaries 

In this manuscript we discuss regularity properties of the boundary 
of the pseudospectrum of a matrix A G C nxn . This boundary turns out 
to be a piecewise smooth curve. Our main concern is how the structure 
of the Schur triangular form of A determines the singular points on this 
curve. 

Let us begin by adapting the results of [2] on general matrix poly- 
nomial to the particular case of the linear monic polynomial. Below, 
Spec A denotes the spectrum of A and || • || denotes the (maximum) 
norm of A as a linear operator in the Euclidean space C n . For a given 
5 > 0, the pseudospectrum of A is the set 

Spe C(5 A:= |J Spec(A + E). 

\\E\\<S 

By construction Spec A = Spec A. For 5 > sufficiently small, 
Spec 5 A consists of "small" connected components around Spec A. As 
5 increases, these components enlarge, collide and eventually intersect 
in various complicated ways and the boundary of Spec^ A, c^Spec^ A, 
becomes more complex. Typically, corners will appear in <9Spec 5 A as 
a consequence of two or more of these components intersecting. As we 
shall see below, the least singular values of -P(A) = (A — A) plays an 
important role in the dynamics of this process. 

Let A = x + iy = (x, y) be in the complex plane (= IR 2 ). The singular 
values of -P(A) are the nonnegative square roots of the n eigenvalue 
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functions of P(X)*P(X). They are denoted by 

si(A) > ... > s„(A) > 0. 

Both the spectrum and the pseudospectra of A are characterised by the 
real- valued function s n : C — ► [0, oo), given by the smallest singular 
value. Indeed, it is well known [9] that for all 5 > 0, 

(1) Spec^ A = {X G C : s n (A) < 5}. 

It will be seen that this characterisation is crucial for the study of 
smoothness properties of the boundary, <9Spec 5 A. They are the subject 
of the next section and the main conclusions are contained in Theorem 
[31 In particular, the notion of fault points (at which the least singular 
value is multiple) is introduced and developed. In Section 3 the Schur 
triangular form is studied and it is shown that, by defining a certain 
equivalence relation on the eigenvalues of A, the general Schur triangle 
can be reduced to a block diagonal form. It is natural to ask whether the 
equivalence classes of eigenvalues generated in this way are independent 
of the the particular unitary transformation used. This is shown to be 
the case when all eigenvalues of A are distinct. 

Section 4 includes a classification and study of singular points on 
the boundary of a pseudospectrum and, finally, Section 5 is devoted to 
a study of these singular points in the case of matrices with size not 
exceeding n — 3. 

2. Regularity properties of the pseudospectral boundary 

Let us first consider regularity properties of s n (A) as a function de- 
fined on the complex plane. Let 

S J :={(x,y, Sj 2 (x,i/))GR 3 : (x,y)eR 2 }. 

Lemma 1. Uj=i ^ s a rea ^ algebraic variety. 

Proof. Define the function 

(2) d(x,y,S) = det [I S - P(x + iy)* P(x + iy)} , (x, y, S) G M 3 . 

As the matrix P(x+iy)*P(x+iy) is hermitian, d(x, y, S) is a polynomial 
in x, y, S with real coefficients, and since 

n 

|J = 2/, SOGR 3 : d(x,y,S) = 0}, 

the result follows. □ 

Real algebraic surfaces usually have strong smoothness properties, 
but self intersections and cusps of various types may also occur as the 
following two examples demonstrate. 
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Example 1. If A = diag [1, e 2in/3 , e" 2i7r/3 ], then 

Spec[P(A)*P(A)] = {|A - 1| 2 , |A - e 2i7r/3 | 2 , |A - e~ 2m/3 \ 2 }. 

Thus Ei U E 2 U E 3 is the union of three paraboloids of revolution with 
minima at 1, e 2l7T ^ 3 and e _2Mr / 3 respectively. These paraboloids intersect 
each other in three different planes parallel to the vertical axis. The 
only point where the three surfaces intersect simultaneously is (0, 0, 1). 

o 

Example 2. If 

/3/4 1 1 \ 
A = 5/4 1 , 
\ -3/4/ 

Ei U S2 U E3 has a conic double point at (0,0,5/16). At this point 
s 3 (0) = s 2 (0) = VE/4. o 

It will be seen below that the occurrence of isolated singularities such 
as that of Example [2] is rare. In this example, the matrix A had to be 
carefully crafted to allow the conic double point around the origin. Any 
slight change in the coefficients of A would eliminate this degeneracy. 

It is well known [4J Theorem S6.3] that, for each j, Sj(j(t)) 2 is a 
real analytic function of t e K whenever 7 : R — > C is analytic. 
Thus, Sj(X) = s*;(A) for j 7^ k on a non-empty open set O only when 
O = C. Therefore two different surfaces Ej can intersect only in a set 
of (topological) dimension at most one (see also [2]). 

To describe these surfaces more precisely, let Pa = P '■ {1> • • • , n} — ► 
{1, . . . , n) (not always onto) satisfy the following properties: 

(a) Ej = Spy). 

(b) Ui=i S i = U"=i S p0) and 

(c) E p(j -) = E p(fc) p(j) =p(/c). 

The map p is a choice of the indexes of the different Ej and only these. 
We denote by thaU) = m (i) = m (p(j)) the (constant) value of the 
algebraic multiplicity of Sj(X) for almost all A G C. 

We are now ready to identify the region of the plane where Sj(X) is 
differentiate. Here and below d$ denotes partial differentiation with 
respect to S. With the determinant function d as in (j2J), let 

G A := {(x,y) G M 2 : d s d l l m ^{x,y,s 2 {x,y)) = 0}. 

Lemma 2. Lei O C R 2 an open set such that On (^U Spec A) = 0. 
TTten s 2 (x,y) is differentiable in O. 
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Proof. This lemma is a consequence of the implicit function theorem 
applied to the function d(x, y, S). Clearly, the polynomial d(x, y, S) is 
continuously differentiable and, since OflSpec A = 0, d(x, y, s?(x, y)) ^ 
for (x, y) G O. These two facts ensure that d l ^ m ^>{x 1 y, S) is also con- 
tinuously differentiable in a suitable M 3 neighbourhood of E 3 - fl(0 X R). 
Since OC\Q A = 0, the hypotheses of the implicit function theorem ap- 
plied to £j- are fulfilled. Since s?(x,y) is the actual "implicit" function 
given locally (and hence globally) by this theorem, the result follows. 

□ 

Note that we can compute Vs? on O by implicit differentiation. 
The sets Q 3 A introduced above are closely related to the set of "fault 
points" introduced in [2J. For j = 1, . . . , n, let 

Fa '■= {( X ,V) '■ s p(j)( x ^y) = s p(k)(x,y) for some k s.t. p(j) ^ p(k)}. 

Here, we call Ta '■= T A ^ ne se ^ °f f au lt points of A. Typically Ta is 
made up of fault lines or curves where Vs n (A) is undefined. 

Let {pi, . . . ,pi} := p({l, • • • , n}). By construction, s Pj ^ s Pk almost 
everywhere and 

d(x,y,S) = (S-s 2 pi r^---(S-sl i r^. 

Also 

dsd y^\x, y ,s) = ^^(s-si j )^- 1 n (s-4 k )^. 

j=l \P % > k^j,l<k<n 

So (x, y) satisfies s Pi (x, y) s Pj (x, y) for all i ^ j if and only if 

d s d 1 / m ^\x,y,s Pi {x,y))^Q. 
It follows from this observation that 

(3) Q\ = T\, j = l,...,n. 

In particular s n (A) will be differentiable outside Tj^. This justifies the 
name chosen for the latter set. 

Clearly [2J, 

(4) <9Spe C(5 AC{zeC : s„(A) = 5}. 

Note that equality in (Tjf]) does not hold in general. In Example 1 above 
the origin is in the right hand set, but it is an interior point of Spe^ A. 

Theorem 3. For any 5 > the boundary of Spec $ A is a piecewise 
smooth portion of an algebraic curve. In particular, it has a finite 
number of singularities. These singularities are either cusps or self- 
intersections. If \q G dSpec$A is a cusp, then A G Ta- 
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Proof. The first part of the theorem is Theorem 7 of [2] . For the latter 
part note that if a singularity occurs at Ao G <9Spec $ A, then either 
A G T a or Vs n (A ) = 0. If Vs n (A ) = 0, then either A is a point of 
self-intersection or Ao G Spec A. Since no eigenvalue of A lies on the 
pseudospectral boundary, the result follows. □ 
Theorem [3] extends to matrix polynomials without much effort. 

3. Refinement of the Schur triangular form 
We denote a Schur factorisation of A by A = USU*, where 

/ail tu ••• ti n \ 

(5) S= ° " 2 : 

t(n— l)n 

\0 ... a n J 

and U is a unitary matrix. Note that Spec A consists of the distinct 
ay, and that S and U are not uniquely defined. The singular values 
of -P(A) are invariant under unitary similarity transformation, so the 
Schur form is an invariant as far as spectrum, pseudospectrum and 
fault points are concerned. 

The characterisation of all possible Schur triangular forms of a gen- 
eral matrix A is certainly beyond the scope of this manuscript. How- 
ever, Theorems [5] and [6] below suggest how this issue can be approached. 
Some preliminary considerations are required. 

Definition 1. Let A = U*SU be a Schur factorisation of A with S as 
in (J5J): 

(a) We write — ctj and say that aj is directly related to otj if one of 
the following conditions holds: 

(i) i = j, (ii) i < j and ty ^ 0, (iii) i > j and tji ^ 0. 

(b) We write ol{ o ctj and say that a» is block equivalent to aj if and 
only if there exists a subset {cxi, . . . , a m } of {1, ... , n} such that 

In this definition we consider that any two eigenvalues on the diago- 
nal of S are "different" even in the presence of multiplicity. The symbol 
"=c=" obviously defines an equivalence relation on the set {ai}™ =1 . 

Let D be the binary matrix obtained from S by preserving the zeros 
of the latter and replacing the non-zero entries by 1. Then D is the 
adjacency matrix (see [3]) of a graph, G. Two diagonal entries of S 
will be block equivalent, o aj, if and only if the nodes i and j of G 
are connected with a path. 



Lemma 4. Let o be any 'permutation of {l,...,n} and S be upper 
triangular, as in (jSJ). Then there is a unitary matrix V such that 



v*sv 



/a<r(i) r i2 
a a {2) 



V o 

and cii o aj if and only if ot. a u) 



Tin \ 



a 



Proof. By writing the permutation a as a product of transpositions 
(cycles of length 2), we see that the case n = 2 can play an important 

«! t 



role. In this case S 



a 2 



with t ^ (so that a% c= a 2 ) and 



there is just one permutation a of interest: (1, 2) — > (2, 1). 

Let a 2 — «i = 6e* e (6 > 0), and a = \t\ and consider the real orthog- 
onal matrix 



W :-- 



a —b 
b a 



Va 2 + b 2 
A little calculation shows that 

a 2 \t\e ie 



w*sw 







a i 



Since t ^ 0, |t|e* e ^ 0, and ct2 =c= oti, as required. 

For the case n > 2 let a k = \tk,k+i\, ®k+i — &k = bke l9k , {b k > 0), and 

let 



1 



bk 



-b k 
a k 



Define 

V k = diag (4-i, W k , In-k-i 
Using a conforming block structure for S, 

Sn SuWk 5*i3 



v;sv k 



o w* k s 22 w k w* k s 23 

533 



=: R 



and 



w;s 22 w k 



ttfe+i t k+ i tk e Wk 
a fc 

By construction, the upper-right entry of the latter matrix is zero if and 
only if tk+i,k = 0. Thus ak — c*k+i (using S) if and only if ak+i — ak 
(using R). 

Since the general permutation a can be expressed as a product of 
"contiguous" cycles as used above, the required transforming matrix V 
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can be expressed as a product of the elementary unitary matrices 14 
introduced above. □ 

We now show that a block diagonal factorisation for A can be ob- 
tained from the relation "=c=". 

Theorem 5. Given any upper triangular matrix S G C nxn there exists 
a unitary matrix U preserving u o " such that 

(6) U*SU = dmg[B 1 ,...,B k ], 

where, for I — 1, ... ,k, 

/5i(0 r 12 (l) ... r lmi (l) \ 

o a 2 (0 '■■ ■ 



\ ... 

and o for all i,j and I. 



r (mi-l)mi (0 



Proof. Since is an equivalence relation, it partitions Spec S into a 
family of equivalence classes. So there is a permutation 



a = {a(l),...,a(n)} 
and there are positive integers l\ < I2 < ■ ■ ■ < h < n such that 
{{a CT (i), • • • , aa(h)}i {a ff (i 1+ i), a CT (« 2 )}, • • • , {a a Q k+1 ), . . . , a a ^}} 

is the partition of Spec S under "=c=". The U found in Lemma H] for this 
a provides the desired conclusion. □ 

In this theorem, <Sj(7) o aj(l) if and only if I = I. The blocks B[ are 
not necessarily unique in the representation U*SU = diag [Bi, . . . , Bk\. 
The diagonal blocks can be permuted, for example. 

Note that the equivalence classes given by "=c=" on the set {ai}2 =1 , 
correspond to the connected components of the graph G introduced 
above. When S is sparse, strategies for computing the block diagonal- 
isation (jH]) may be based on this observation. 

It is natural to ask whether Definition [T] is independent of U and S 
in the Schur decomposition for A. We obtain a positive answer if all 
eigenvalues of A are simple. 

Theorem 6. Suppose that A has no multiple eigenvalues, i.e. aij 7^ aj 
for i 7^ j . Then the equivalence classes determined by the equivalence 
relation "o " of Definition [I] are independent of the Schur factorisation 
chosen for A. 
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Proof. This result is a direct consequence of the following observation. 
If S and T are two triangular matrices with the same diagonal entries 
Su = ta = ccj and with aj ^ otj for i / j, such that U*TU = S 
for a suitable unitary matrix [/, then U must be diagonal. See [8j 
Theorem 2.3] or [7]. □ 

The results reviewed at length in [8] suggest that extension of this 
theorem to admit multiple eigenvalues of A would be difficult. 

4. The singular points on the pseudospectral boundary 

As we shall see next, the block diagonalisation found in Theorem [5] 
provides a natural classification of the singular points on the boundary 
of the pseudospectra of A. 

Let S be a Schur triangular form of A and B\ be the blocks corre- 
sponding to the diagonalisation of S given in ([6]). It is clear from the 
definition that the singular values of P(X) = A — A are those of (X — Bi) 
for I = 1, . . . , k. Then it follows from ([TJ that 

k 

(7) Spec 5 A = IJ S-pec s Bi, for all 5 > 0. 

1=1 

This decomposition motivates the following classification of singular 
points on the boundary of the pseudospectrum: 

Definition 2. We say that Ao G dSpec s A is a 

- stationary point if A G" Ta but Vs n (A ) = 0, 

- essential fault point if Ao G Tbi for some / = 1, . . . , k, 

- regular fault point if Ao G ^4 but it is not an essential fault 
points. 

Clearly a stationary point of A must also be a stationary point of B\ 
for some (not necessary unique) I = 1, . . . , k. 
If A G Ta, then either A G [jf =1 ^b, or 

it 

Ao G (<9Spec 5J B fe n dSpec^) \{jT Bl 

1=1 

for some k ^ I. In the former case, Ao is an essential fault point and in 
the latter it is a regular fault point. 

Since it is formed as a consequence of two pseudospectra of different 
blocks Bi intersecting, whenever non-empty, the regular portion of Ta 
is expected to be of topological dimension 1. On the other hand, we 
shall argue in Section [5] that blocks of small size have only a limited 
number of essential fault points. 



The role played by the classification just introduced in the dynam- 
ics of the pseudospectrum of A can be better visualised by means of 
concrete examples. However, let us first establish two elementary con- 
sequences of (J7J). 

Corollary 7. If A is a normal matrix, then all singular points on 
dSpec s A are regular fault points. Furthermore, Ta is the Voronoi dia- 
gram associated to Spec A. 

Proof. The proof is straightforward. See Example [1] for illustration. 

□ 

Corollary 8. Let A be unitarily similar to a bi-diagonal matrix. Then 
Ta = if and only if no entry in the off-diagonal of A vanishes. 

Proof. Use the fact that (A — B)*(X — B) is a tri-diagonal symmetric 
matrix. See [TUJ §5.36]. □ 

Corollary [H] implies that bi-diagonal matrices have no essential fault 
points. 

Example 3. A detailed analysis is made of a problem with two diag- 
onal blocks (k — 2); see also Figure [2j Let 

/3 N 
A= -1 1 

V° \ 

Any pseudospectrum of A is the union of those of the blocks 

-1 f 



B\ = 3 and B2 



1 



The (only) singular value of [x + iy — B\ ) is \J {x — 3) 2 + y 2 , so Spec^-Bi 
is a disk centred at (x, y) = (3, 0) of radius 5. The singular values of 
B 2 are 

(j \ V 2 
3 2 2 ^ v/5 + 20x 2 + V X 
~ + x +y ± 

It is straightforward to see that the least singular value s_ (x + iy) = 5 
if and only if 

(x 2 + y 2 ) 2 - 2(1 + 5 2 )x 2 + 2(1 - 5 2 )y 2 + 1 - 36 2 + 5 4 = 0. 

This shows that Spec s B 2 is a spiric section for all 5 > 0. It is also 
straightforward to see that Ta is a portion of the hyperbola given 
explicitly by 

T A = {{x,y) : 31 x 2 - y 2 - 90 x + 55 = 0, x > 1}. 
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Any singularity on dSpec^ A which arises as a consequence of <9Spec 5 B\ 
intersecting with dSpec^f^, will be a regular fault point. Also, the 



curve 9Spec 5 B 2 has a stationary point at the origin when 5 = 




Note that a further self-intersection occurs on dSpec^ A at (x, y) = 
(«±§3£ o) for 5 = since the boundaries of Spec s B x and Spe C(5 B 2 

touch at this point. The pseudospectrum of A will consist of three con- 
nected components for < 5 < \J , two components for y^ 3 "^ < 
5 < 48 ~^ and a single component for all 5 > 48 ~^ ■ o 

The next observation may be relevant in the effective design of a 
corrector step in path-following algorithms for tracking <9Spec 5 A. The 
second part asserts that the corners in SSpec^ A at regular fault points 
will always be re-entrant. 

Lemma 9. Let Xq be a regular fault point in dSpec s A. Assume that 
the curve dSpec s A fails to have a tangent line at Xq and does not self- 
intersect at this point. Then we can always find an a > it, depending 
only on Xq, satisfying the following property: If0</3<a, there exists 
r > and < 7 < 2tt such that the sector 

{X = X + pe i(9+7) : 0<6<f3, < p < r} C Spec 5 A. 

Proof. This is a consequence of the fact that a regular fault point of 
<9Spec 5 A is formed by two different intersecting connected components. 

□ 

Theorem 10. Ifoti, ctj G Spec A are in the same connected component 
of C \ Ta, then a, o atj . 

Proof. It suffices to show that if oti <f? ctj, then any continuous trajec- 
tory : [0, 1] — > C such that 0(0) = a« and 0(1) = ay, intersects 
J~s ( see ©)• We achieve this by noticing that a-j G diag[5/] and 
ctj G diag[5 m ] for I 7^ m, and applying the mean value theorem in- 
ductively to pi — pj , where 

p k (t) := min.sing.val(0(t) - B p ), a k G diag [B p ]. 

□ 

The converse of Theorem [TU] does not hold in general. It is easy to 
construct examples where «j o ctj, but a-j, otj both belong to different 
components of C \ Ta- One such example is the following. 
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Example 4. If 



A = 



( 


-1 


i 





\ 







i 
















—i 


1 


V 











* / 



T A = {A G C : Re A = ±Im A}. o 



5. Essential fault points and the Schur structure of 

small matrices 

Determining the structure of J- a for a given matrix A, is typically 
involved. As we confirmed in the previous section, the set of fault 
points can be empty or consist of a single point, but it can also be a 
complicated set such as a Voronoi diagram. By virtue of (J7J), the set of 
regular fault points is completely characterised once the pseudospectra 
of each of the blocks in the diagonal factorisation ([6]) are known. In 
this section we show that, for sufficiently small matrices, the number 
of essential fault points is finite. 

It is easy to characterise the set of fault points of a 2 x 2 triangular 
matrix 

"-(;;)■ 

If a =c= /3, then Tb = 0- If o. ft (3, then 

T B = {A G C : |A-a| = \\-/3\} 

for a f3 and Tb = for a = (3. By virtue of Corollary [HI no fault 
point of a 2 x 2 matrix can be essential. 

The 3x3 case is more involved. This is illustrated in the following 
example (which is a generalisation of Example [2]) . 

Example 5. See Figure [31 Let 

(a 1 l\ 

A = 5/4 1 , 
\0 c) 

where \a\ = \c\ = 3/4. Then Ta — {0} and a o 5/4 o c. Whenever 
G dSpec s A, it will be an essential fault point. This can only occur 
at 5 = s 3 (0) = y/E/4. 

The complicated dynamic of the pseudospectral boundary as we 
move the parameters a and c on the circle with radius 3/4, is illus- 
trated in Figure [3j There we show the evolution of the essential sin- 
gularity when a = c = 3e ie /4 for 6 = kir/A, k = 0,1,2,3,4. When 
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k = 0, Spec ^4 A is connected and there is clear indication of a re- 
entrant corner on the left side of the boundary. The re-entrant angle 
at this corner becomes more pronounced when k = 1. For k = 2, 3 
the concavity prevails. However there exist critical Q\ G (7r/4, 7r/2) and 
9 2 G (37r/4, 7r), such that G" <9Spec^5/ 4 v4 for 6 G (0i,#2)> and hence 
there is no corner on dSpec^^A At the final stage k — 4, the pseu- 
dospectrum now consists of two connected components, the concavity 
on the left side has vanished, however a new re-entrant corner forms in 
a different part of the boundary. o 



More generally, let 



(8) 



B 




When atj ^ for some j ^ k, B can be reduced to a block diagonal 
form of smaller size, so the set of fault points is characterised by using 
(JZj). If ax o ot2 =c= «3, every point in Tb is an essential fault point and 
we have the following result. 

Theorem 11. Let B be as in (JSJ) and assume that no two of the upper 
triangular entries r,s,t vanish simultaneously (i.e. ati o a 2 o a 3y ). If 
rst = 0, then Tb = 0- Ifrst ^ 0, then Tb is either empty or consists 
of a single point. In the latter case, 



(9) 



J~b C {A G C : arg(a 2 — A) = arg(rst)}. 



Proof. If rst = 0, the claimed assertion is a consequence of Corollary [8] 
and a suitable permutation of the rows and columns of B. 

(a± — A), b = (a 2 — A), c = (a 3 — A), and assume that rst ^ 0. 



Let a 
Let 



so that 



and 



u 



and 



B- A 



(B - X)*(B - A) 



Bx 


u 





c 



B{B t 


V \ 


V* 


\\u 2 \\ + 


c 


•2 1 



where v = B*u. By the Cauchy interlacing theorem, if B*B has a 



double eigenvalue a 
of B{Bx. 



Sjj then cr is also the minimal eigenvalue 
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Since r ^ 0, B\Bi has only simple eigenvalues. Let 7^ e G C 2 be an 
eigenvector such that B\B\e = ere. A straightforward argument shows 
that, if a is a double eigenvalue of B*B, then »le. 

N ° W ' W= (r/+5t)' 

2a = |a| 2 + |6| 2 + |rf - v^M 2 + N 2 + M 2 ) 2 - 4|a| 2 |6| 2 and 
_ f\b\ 2 -\a\ 2 + \r\ 2 - y/(\a\ 2 +\b\ 2 +\r\ 2 ) 2 -4\a\ 2 \b\ 2 \ 

6 v ~ 2a¥ J ' 

Note that all these three quantities depend on A. If v*e = 0, then 

q{\) = s(|6| 2 -|a| 2 -|r| 2 - v/(|a| 2 + |6| 2 + |r| 2 ) 2 - 4|a| 2 |6| 2 ) - 2rt6 = 0. 

As the coefficient of s in the above expression is real, (02) is guaranteed. 

Let 7(t) = a 2 + te i£LTg ( rst \ — 00 < t < 00, be a parameterisation of 
the line where Ta lies. Then q("f(t)) = qi(t) — \J q 2 (t) where qi depends 
linearly in t and q 2 is a quadratic polynomial in t. Moreover, {q\{t)) 2 
and q-i (t) have the same coefficient of order 2 in t. Thus (?(7(t)) can 
only vanish at no more than one t- value, t\. Since 7(^1) is the only 
possible essential fault point of B, the proof is complete. □ 

In Example [5] the line described by the right hand expression in (J9j) 
is the real axis and Ta is the origin. 

It is natural to expect that the argument presented in the proof of 
Theorem [11] can be extended inductively to blocks of larger size. We 
have not explored this possibility in much detail. However our observa- 
tions lead us to conjecture that, perhaps, the number of essential fault 
points is always finite for matrices of any size. This issue certainly 
requires further investigation. 
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FIGURE 2. Example [3l The thick lines are dSpec^A for 5 — 
2/5, s/*=^, 48 ~^ ,8/5. The thin line is T A . 



Figure 3. Evolution of the essential singularity on the pseu- 
dospectral boundary for the family of matrices given in Example[5l 
Here a = c = 3e ife7r / 4 /4 for k = 0, 1, 2, 3, 4 (top to bottom). On the 
left side we depict SSpec^A for 8 = On the right side 

we show details of the pseudospectral boundaries near the origin 
for S close to the critical value V5/4. 



